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m , ABSTRACT 

The massive clump GlO.6-0.4 is an OB cluster forming region, in which multiple UC 
Hll regions have been identified. In the present study, we report arcsecond resolution 
observations of the CS (1-0) transition, the NH3 (3,3) main hyperfine inversion transi- 
tion, the CH3OH J=5 transitions, and the centimeter free-free continuum emissions in 
^ this region. The comparisons of the molecular line emissions with the free-free contin- 

^ • uum emissions reveal a 0.5 pc scale massive molecular envelope which is being partially 

^ dispersed by the dynamically-expanding bipolar ionized cavity. The massive envelope 

is rotationally flattened and has an enhanced molecular density in the mid-plane. In 
the center of this massive clump lies a compact (<0.1 pc) hot (^100 K) toroid, in which 
a cluster of 0-type stars has formed. 

This overall geometry is analogous to the standard core collapse picture in the low- 
mass star forming region, with a central (proto-) stellar object, a disk, an envelope, and 
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a bipolar outflow and outflow cavity. However, GlO.6-0.4 has a much larger physical 
size scale (<0.1 pc for typical low-mass star forming core). Based on the observations, 
we propose a schematic picture of the OB cluster forming region, which incorporates 
the various physical mechanisms. This model will be tested with the observations of 
other embedded OB clusters, and with numerical simulations. 



Subject headings: stars: formation — ISM: evolution — ISM: individual (GlO.6-0.4) 

1. Introduction 

A M>10 Mq massive star has a short Kelvin-Helmholtz (cooling) timescale, and starts the 
nuclear burning while it is still deeply embedded in the dense molecular clump. During this earliest 
evolutionary stage, the massive stars can continue to increase their stellar mass via accretion of more 
dense molecular gas (Kahn 1974; Garay and Lizano 1999). Meanwhile, these massive stars produce 
large amounts of strong ionizing photons, and create the Ultracompact (UC) Hii Regions. The 
UV radiation (Krumholz, Mckee, Klein 2005; Peters et al. 2010), the pressure of the ionized gas 
(Keto 2002b; Krumholz et al. 2009; Peters et al. 2010) , and the strong stellar wind directly interact 
with the surrounding molecular gas. In addition, as massive stars typically form in clusters, their 
environments may be also influenced by their interactions both radiatively and dynamically. These 
interactions potentially can disperse the accretion flow, unless the system is confined or trapped by 
the high density molecular gas with the right physical profiles (example solutions see Keto 2002b, 
2003, 2007; also see Galvan-Madrid et al. 2008 for an observational evidence of accretion). This is 
the greatest contrast to the formation of the low-mass stars, where accretion only occurs before the 
pre-main sequence phase when the surrounding gas is neutral. To understand how massive stars 
finally attain a high stellar mass (e.g. >20 Mq), it is necessary to understand the accretion after 
the formation of an UC Hll. One key factor in determining the relative importance between these 
feedback mechanisms as well as gravity, and hence the subsequent evolution of the dense molecular 
core/envelope, is the overall geometry /morphology of the entire system. For example, the density 
distribution of the molecular gas around massive stars determines how the molecular structures are 
self-shielded from the ionizing photons; and the overall geometry of the system determines how the 
radiation and the ionized photons can leak out of the system. 

Owing to their typical large distances (a few kpc), the molecular structures around the UC Hii 
regions, are not well resolved in previous observations. How the accretion fiows around individual 
stars and in the entire cluster is unknown. Hence how the OB stars gain mass via accretion after 
nuclear burning begins, is also unknown. We perform new interferometric observations toward a 
well-studied O-type cluster forming region GlO.6-0.4, in order to resolve the molecular structures. 
The high column density and high brightness temperature of this source allow us to follow the 
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detailed structures of the molecular gas over an extended region, at arcsecond resolutions. Here we 
emphasize that the subject under consideration is the molecular and the ionized gas flow around a 
cluster of OB stars, rather than the accretion around a single massive (proto-) stellar object. This 
is due to the still insufficient spatial resolution. 

The UC Hii region complex GlO.6-0.4 is at a distance of 6 kpc (Caswell et al. 1975; Downes 
et al. 1980), and is still deeply embedded in a dense molecular core/envelope. The brightest UC 
Hii region is extremely luminous (integrated flux is 2.6 Jy at 23 GHz) (Ho & Haschick 1981; Sollins 
et al. 2005; Sollins & Ho 2005), and contains an OB cluster within a 0.05 pc region. Detections 
of water and OH masers (Genzel & Downes 1977; Ho & Haschick 1981; Ho et al. 1983; Hofner & 
Churchwell 1996; Fish et al. 2005 ) indicate that the formation of massive stars is still ongoing. 
The early NH3, CS and C^^O line observations (Ho & Haschick 1986; Keto, Ho, k Haschick 1987; 
Keto, Ho, & Haschick 1988; Omodaka et al. 1992; Ho, Terebey, & Turner 1994 ) unveiled a 0.5 pc 
scale molecular envelope. These observations also suggested that the dynamics of the majority of 
the molecular gas have not yet been strongly disturbed by the YSOs or stellar activities. At the 
0.5 pc scale, the general motion of the molecular gas appears to be dominated by gravity, showing 
rotation along a flattened geometry (Liu et al. 2010). The inner part at the 0.1 pc scale appears 
to rotate faster while infalling towards the UC Hii region (Keto, Ho, &: Haschick 1987; Keto, Ho, 
& Haschick 1988; Liu et al. 2010; Beltran et al. 2011). Recombination line studies further suggest 
that the accretion flow continues across the Hll boundary, while evidence of outflow can be seen in 
the motions and structures of the Hii region at the sub-arcsecond angular scale (Keto 2002; Keto 
& Wood 2006). Structures found in projection against the UC Hll continuum emission from NH3 
optical depth studies, suggest accretion via a very clumpy and geometrically thick rotating core at 
the 0.1 pc scale (Sollins & Ho 2005). 

In this paper, we present the new arcsecond resolution (1" corresponds to 0.03 pc at the 
distance of GlO.6-0.4) molecular line (CH3OH, NH3, CS) data, together with very sensitive and 
0".5 resolution 3.6 centimeter continuum observations. The linear scale of the synthesized beams 
is smaller than the natural scale for self-gravity (Jeans length ~0.1 pc for temperature ~20-30 K 
and H2 number density ~10^ cm~^). We see an unprecedented detailed morphology of the massive 
accretion flow which clarifies the geometrical relation between the hot core and the envelope. With 
the aid of the new results and the discussions on this target source in previous publications, a more 
thorough understanding is achieved. 

In massive star forming regions, the NH3 inversion transitions are usually regarded as reliable 
molecular tracers for 20-100 K gas. The main hyperfine components of these NH3 lines are optically 
thick and may mask the detailed density structures. The satellite hyperfine components are more 
optically thin (see Ho & Townes 1983 for a review). However, in most cases, the satellite hyperfine 
components of the NH3 (J, K) = (1,1) inversion transitions are blended with the main hyperfine 
component; the satellite hyperfine components of the higher (J,K) level transitions are much weaker 
and are difficult to be observed with high resolutions. With a factor of ~100 higher Einstein A 
coefficient, the selected CH3OH transitions potentially allow us to probe the denser regions. We si- 
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multaneously observed multiple CH3OH transitions which have comparable Einstein A-coefficients 
(and therefore critical densities), but trace a broad range of upper-level-energy (E„p). These ob- 
servations allow us to distinguish the hottest molecular gas directly associated with the UC Hii 
region from the cooler envelope. 

The low E„p and the modest critical density of the CS (1-0) transition provide sufficient optical 
depth, which allows us to trace the molecular gas to a much more extended region. However, in the 
case of GlO.6-0.4, significant local structures can still be identified without being severely affected 
by the sidelobes of the dirty beam and the missing fiux. The high optical depth of CS (1-0) 
also provides better signal-to-noise ratio to the geometrically thinner expansional signatures, and 
provides information for the on-going feedback processes. 

The observational results are presented in Section [3l In Section HJ we discuss the physical 
implication of the overall geometrical picture of the region. In addition, we estimate the feedback 
from various physical processes (radiative pressure, ionized gas pressure, stellar wind) and compare 
with our observational measurements. A brief summary of the results are provided in Section [5l 



2. Observations and Data Reduction 

2.1. The Molecular Line Observations 

We observed the CH3OH J=5 transitions in the 230 GHz band using the Submillimeter Array 
(SMA)E|; we observed the NH3 (3,3) hyperfine inversion transitions in the K-band using the NRAOg 
(Expanded) Very Large Array (VLA/EVLA) in the C-configuration; and we observed the CS (1- 
0) transition using the VLA/EVLA in the DnC-configuration. For more on the SMA and its 
specifications, see Ho, Moran &; Lo (2004). The size of the synthesized beams in these observations 
are comparable. Continuum emissions are averaged from the line-free channels and then subtracted 
from the line data. The properties of the selected molecular transitions are listed in Table 1; the 
instrumental parameters are summarized in Table 2; the observational settings are provided in 
Table 3. The CH3OH 5(-2,4)-4(-2,3) E and CH3OH 5 (2,3)-4(2,2) E lines are closely blended. The 
frequency separation of these two lines is only marginally resolved. Considering that these two 
lines trace similar physical conditions, and the line-widths are broad in the observed region, we 
regard the blended lines as one single line in this research, and denote it henceforth as CH3OH 
5(2,3)-4(2,2). 



^The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the Academia 
Sinica Institute of Astronomy and Astrophysics, and is funded by the Smithsonian Institution and the Academia 
Sinica. 

^The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under 
cooperative agreement by Associated Universities, Inc. 
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Basic calibrations are carried out in the MIR package for the SMA data; and are carried out 
in the AIPS package for the the VLA/EVLA data. All data are additionally self-calibrated in the 
AIPS package. Imaging is also carried out in the AIPS package. 

2.2. The Centimeter Continuum Observations 

We observed the X band continuum emission toward GlO.6-0.4 in the VLA A-configuration 
including the VLB A Pie-Town antenna on 2005 January 2, and in the VLA/EVLA C-configuration 
on 2009 July 27. We observed 1331+305, 0319+415, and 1820-254 as absolute flux, passband, and 
gain calibrators. 

The basic calibrations, self-calibration, and imaging of these data are carried out in the AIPS 
package. We combine the A-array+Pie-Town visibility data with the C-array visibility data, 
yielding a 0".5x0".4 synthesized beam with a position angle of -43.2°. The observed RMS noise 
of the 3.6 cm continuum image is about 0.24 mJy/beam (~21 K in brightness temperature), and 
the peak brightness temperature is about 8700 K. We note that the 3.6 cm continuum emission is 
free-free emission based on spectral index measurements. Hence its high brightness temperature 
indicates that the emission is moderately optically thick. 

3. Results 

The observing results are separately presented in two sections. In Section 13.11 we focus on 
discussions of the warm and dense molecular gas, mainly traced by the NII3 and the CH3OH 
emissions. In Section 13. 2( we examine the dynamics of the more extended components, based on 
the observations of the CS (1-0) transition. Cross comparisons of the molecular lines have been 
omitted, when there are no additional significant physical implications. 

3.1. The NH3 and the CH3OH Emissions 

3.1.1. The Distribution and the Dynamics of the Dense Gas 

Figure [1] shows the velocity integrated maps of the CII3OH 5(0,5)-4(0,4) A+ transition, the 
blended CH3OH 5(-2,4)-4(-2,3) E and CH3OH 5(2,3)-4(2,2) E transitions, and the CH3OH 5(- 
3,3)-4(-3,2) E transition. In Figure [21 we present these three images in RGB colors, and an overlay 
with the velocity integrated map of the NII3 (3,3) main hyperfine line. In the projected area of 
the UC Hii region, the NH3 (3,3) main hyperfine line is detected in strong absorption (see Keto, 
Ho k Haschick 1987; Keto, Ho, & Haschick 1988; Solhns et al. 2005; Liu et al. 2010 for more 
descriptions) and is presented with much larger contour intervals. We see that: 
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• The highest excitation CH3OH transition (blue) exclusively traces a compact (~0.1 pc) flat- 
tened hot toroid in the center (~300K, Klaassen et al. 2009; ~87 K, Beltran et al. 2011). The 
lower excitation CH3OH transitions trace the warm (~50 K; Keto, Ho &: Haschick 1987) and 
dense molecular gas to an extended (~0.5 pc) region. The ratio of the velocity integrated flux 
of the CH3OH 5(0,5)-4(0,4) A+ transition to that of the CH3OH 5(-3,3)-4(-3,2) E transitions 
in the hot toroid region is about 0.56. In the optically thin limit, this line ratio reflects an 
averaged excitation temperature of ~500 K. If the lower excitation line is saturated owing to 
high optical depth, the excitation temperature would be overestimated. 

• We see the complicated structures from the extended CII3OH emissions. The majority of the 
CH3OH emission has a southeast-northwest distribution. Two low column density cavities 
are clearly seen northeast and southwest of the hot toroid. 

We note that the size scale of the region where we detect significant emissions is small as 
compared with the angular size scale of the SMA primary beam (~l') and that of the VLA primary 
beam (~2'). Since the detected structures are well centered within the primary beams, they are not 
affected by the primary beam attenuation. We therefore do not perform primary beam corrections 
to the maps. The shortest baselines of the CII3OH and the NII3 data sets (Table [2]) correspond to 
the maximum angular size scale of ~35" and ~90", respectively. Hence the absence of structures 
with angular size scales >20" suggests that the compact appearance of the emissions is physical in 
nature, and is not caused by sampling defects. 

From Figure we see that most structures traced by the CH3OH emissions have their NH3 
counterparts; however, the NH3 emission shows more significantly detected structures than the 
CH3OH emissions. A similar effect appears in the velocity domain. Figure [3] shows the position- 
velocity (PV) diagrams of the NH3 (3,3) main hyperfine line and all selected CH3OH J=5 transi- 
tions. The PV cut is centered on the brightest free-free continuum soured at the coordinates of 
R.A. = 18''10'"28.64^ and Decl = -19°55'49.22" with a position angle PA = 140° (the positive angle 
is defined in the usual counterclockwise fashion). From these PV diagrams, we see that: 

• The CH3OH 5(-3,3)-4(-3,2) E transition has the highest upper-level energy E„p of 96.9 K. 
The PV diagram of this transition reveals the fast rotating hot toroid. 

• The CH3OH transitions with lower upper-level energies trace the molecular gas in a more 
extended region. A significant velocity gradient is also seen in the extended emission. 

• From the left most panel of Figure [3l we see that the dominant emission of the NH3 (3,3) main 
hyperfine inversion line shows a good agreement with the CH3OH emissions. The NH3 (3,3) 
absorption line consistently traces the fast rotating hot toroid as seen in the PV diagram of the 
CH3OH 5(-3,3)-4(-3,2) E transition. Additionally, the NH3 (3,3) emission shows significant 



■^This is exactly the same PV cut which shows the maximum velocity gradient (Liu et al. 2010). 
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detections of broad line emission features (e.g. at 7" and 11") and some diffuse broad line 
emission (e.g. around -3"- -11"). 

The NH3 molecule has a stable abundance over a large range of densities and temperature 
without depletion onto grains. Hence its emission lines can be used to trace the dynamics in the 
bulk of the dense molecular accretion flow. The detected velocity gradients can be explained by the 
marginally centrifugally supported motion of the molecular gas, in response to the gravity of the 
enclosed molecular and the stellar mass (see Liu et al. 2010 for more discussions). The projected 
area and the rotational motions traced by the CH3OH are highly consistent with those which are 
traced by the NH3 emissions. This suggests that these two molecules are tracing dense gas in the 
same region. Nevertheless, there are differences in their detailed brightness distributions which sug- 
gest that in this particular source, the CH3OH emissions are tracing different excitation/chemical 
conditions (e.g. temperature), in the same flow traced by NH3 (see also the discussions for the 
CH3OH emissions in Liu, Ho & Zhang 2010). In addition, the abundance of the CH3OH molecule 
may also be affected by outflow. 

The NH3 (3,3) main hyperfine inversion transition has comparable upper-level energy (E^p/k 
= 124.5 K) with the highest excitation CH3OH 5(-3,3)-4(-3,2) E transition (E^,p/k = 96.9 K), 
but in contrast it shows an extended distribution. One explanation is that the NH3 molecule has 
a higher abundance ([NH3]/[H2] = lO"'^; Keto, Ho & Haschick 1987) than the CH3OH molecule 
([CH30H]/[H2] = 2-10"^; Takakuwa, Ohashi, & Hirano 2003), and thus becomes optically thick and 
shows strong emission. The excellent match between the NH3 emission and the lowest excitation 
CH3OH emission suggests that the emissions of both molecules trace the real physical structures. 
We know that for the dense gas with temperature around 50 K, the statistical weight compensates 
the Boltzmann factor, and therefore the population of NH3 at the (3,3) level can be comparable to 
the population at the lowest excitation (1,1) level. The CH3OH emissions and the NH3 (3,3) main 
hyperfine emission may both trace a population of molecular gas which is significantly warmer and 
denser than the ambient molecular gas. We call this warm (>30 K) and dense (>1(P cm~^) region 
the massive envelope hereafter. The previous observations on C^^O (2-1) unveiled a parsec scale 
filamentary structure in the south of the massive envelope (Ho, Terebey, &: Turner 1994). The 
distribution of the molecular mass on an even larger scale has to be further examined by sensitive 
single-dish maps. 

The selected CH3OH transitions have higher Einstein A-Coefhcients than the NH3 (3,3) main 
hyperfine inversion transition, and potentially only trace the denser components. The observations 
in the clearest example of photon-dominated regions (PDR), the Orion Bar, further suggest that 
the CH3OH emissions trace the clumpy materials, which have an H2 volume density higher than 
10^ cm~^, and are self-shielded from photo ionization (Leurini et al. 2010). The abundance of 
the CH3OH molecule can be greatly enhanced if the molecular gas is shocked by high velocity 
outflows with velocities greater than ~20 kms~^ (see discussions in Takakuwa, Ohashi, & Hirano 
2003). We therefore cannot rule out the possibilities that zones with locally enhanced CH3OH 
abundance/emission are distributed over the entire region. 
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Part of those extended NH3 emissions can also be attributed to the externahy heated lower 
density surfaces of the extended structures. Other possibilities will have to be investigated by higher 
resolution observations which define the detailed temperature and density profiles of individual local 
structures. 

3.1.2. Comparison with the Distribution of the Free-Free Continuum Emissions 

In Figure m we compare the distribution of the molecular gas traced by the CH3OH transitions 
with the distribution of the ionized gas traced by the 3.6 cm free-free continuum emission. By 
comparing the left panel of Figure [2] with this figure, we find that the 3.6 cm free-free continuum 
emissions peak at the hot toroid. The distribution of the extended free-free continuum emission 
fits into the projected area of the low column density cavities seen in the CH3OH RGB map. 

These observed features suggest an overall geometry of a core-envelope-cavity system. This 
overall geometry resembles the standard disk-envelope system in the low mass (<1 Mq) star forming 
region with bipolar outflow and outflow cavity, except on a larger and more massive scale. This 
overall geometry remains even after the embedded OB cluster starts the nuclear burning, and other 
physical processes start to affect the dynamical evolution (Section [T|). We emphasize that the 
subject under consideration is the formation and accretion of a cluster of OB star. By comparing 
the flux of the ionizing photons and the bolometric luminosity. Ho and Haschick (1981) suggested 
that a few massive stars (from 06.5 to BO) have already been formed and are embedded in the 
central 0.1 pc region (see also the discussions in Keto 2002, Sollins &: Ho 2005, and Keto & Wood 
2006). This cluster of OB stars contains ~200 M0 of stellar mass in total. As pointed out by 
Sollins and Ho (2005), the source GlO.6-0.4 stands distinct from the cases where a single massive 
star dominates at the center (e.g. IRAS 20126, IRAS 18089, 0192.16, etc. Zhang et al. 1998; 
Beuther et al. 2004; Shepherd et al. 2001), which might have formed in a process similar to 
low-mass stars (Keto &: Zhang 2010). 

In GlO.6-0.4, a parsec scale massive envelope (>2500 Mq over the central 30" region; Ho, 
Terebey & Turner 1994) appears to be undergoing a coherent motion in response to the gravity 
of the enclosed mass and the specific angular momentum. In this case, the mass of the embedded 
OB stars is a less significant fraction (<10%) as compared to the molecular mass in the envelope. 
The fact that a concentration of OB stars formed in the small central 0.1 pc region suggests that 
the global contraction of the massive envelope may have been efficient. We note that in another 
OB cluster forming region G20.08-0.14N (L~6.6T0^ L0, d~12.3 kpc), the emissions of the NH3 
molecules and various hot core tracers also show a similar geometry and kinematic profile with 
those in GlO.6-0.4 (Galvan-Madrid et al. 2009). 
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3.2. The CS (10) Emissions 

3.2.1. Distribution 

Figures [5] and [6] show the velocity integrated (moment 0) map of the CS (1-0) emission, the 1.3 
mm continuum image, and the 3.6 cm continuum image. In the 3.6 cm continuum image (Figure 
[6]), we see three isolated UC Hii regions (from east to west, UC Hii region A, B, C hereafter), of 
which we also detect their free- free continuum emission at 1.3 mm (Figure O right panel). 

The visibilities of the CS (1-0) data cover a range of uv-distances of 4-245 kX, which cor- 
respond to ~52"-0".8 in angular scale. The lo" rms in the CS velocity integrated map is ~77 
mJy/beam-kms~^ (24 K-kms~^). The achieved brightness sensitivity is comparable to the theo- 
retical rms noise (also achieved in the channel maps). With the presence of the bright UC Hii 
regions, the detection of CS (1-0) is locally biased by the continuum subtraction. Structures dis- 
tributed in front of the UC Hii regions will be detected in emission line only if they emit brighter 
CS (1-0) emission than the local free- free continuum background; otherwise, they will be detected 
in absorption line. We omit plotting the negative flux due to absorption to enhance the contrast 
of emission features in the gray scale. Structures behind the UC Hii region are not affected by 
the continuum subtraction and are always detected in emission line, unless the 49 GHz continuum 
emission is optically thick. In regions where the brightness temperature of the free-free continuum 
emission is comparable with the excitation temperature of the molecular gas, the absolute values 
of the flux density of the emission line and the absorption line are comparable. Smoothed by the 
finite synthesized beam, the emission line and the absorption line signatures may cancel each other, 
and reduce the signal-to-noise ratio. 

From Figures [5] and [6] we see: 

• The brightest regions of the CS (1-0) emission agrees well with the 1.3 mm continuum emis- 
sion. (The significant deficit of CS (1-0) emission in the middle of the maps is due to the 
subtraction of the bright free- free continuum in the UC Hii region.) 

• We see an extended dense flattened structure with a size scale of ~0.5 pc. The orientation is 
indicated by the red dashed line, which has a position angle of 135°, and is consistent with 
the orientation of the central 0.1 pc hot toroid (Liu et al. 2010). 

• Northeastern to the dense flattened structure, the extended CS emission shows filamentary 
structures. The size scale of those filamentary structures is a few times of 0.1 pc. Southwestern 
to the dense flattened structure, we see an upside-down V shaped cavity-wall signature. 

The 1.3 mm continuum emission arises predominantly from the free- free continuum emission 
in the central 2"-3" UC Hii region; in the extended region, the 1.3 mm continuum emission traces 
the thermal dust emission (Liu et al. 2010; Liu, Ho & Zhang 2010). The thermal dust emission is 
usually regarded as a reliable tracer of the molecular hydrogen. The general agreement between the 
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CS (1-0) image and the 1.3 mm continuum image indicates that CS (1-0) also traces the molecular 
hydrogen. However, in the optically thick case, the velocity integrated flux of CS (1-0) can be biased 
by saturation and/or foreground absorption, and also by the non-uniform excitation temperature. 
These explain some inconsistency between the CS (1-0) velocity integrated image and the 1.3 mm 
continuum image. By inspecting the position-velocity diagrams (Section (322]), we do not see a gap 
of emission line, and therefore we argue that the self-absorption does not severely affect the overall 
picture of geometry. The excitation temperature only changes the relative strength of emissions, 
and should not greatly alter the overall picture either, as long as the excitation does not change by 
a very large factor. 

The distribution of the CS (1-0) emissions is consistent with the overall geometry suggested 
by our observations of the CH3OH J=5 transitions. However, instead of seeing a clean biconical 
feature as in the CH3OH maps, we detect significant CS (1-0) emission in the projected area of 
the biconical cavity. This can be due to the fact that the CS (1-0) transition has a much lower 
upper-level-energy (Eup = 2.35 K), and a factor of ~50 lower Einstein A-Coefficient. These two 
factors make CS easily excited in the outer layer of the cavity-wall with lower temperature (and 
molecular gas density) . This property allows us to robustly detect the dynamics of the cavity-wall, 
and explicitly measure the expansional velocity. 

Our CS (1-0) maps additionally unveil an unprecedentedly detailed morphology in the massive 
cluster forming clump. The velocity integrated map apparently shows abundant structures with 
a size scale of a few arcseconds (0.03 pc). In figures El we present an example channel image at 
visr= -7.8 kms~^. In this channel, we do not detect compact (2"-3") absorption features towards 
the brightest UC Hii region A, which alleviates the potential confusion from the sidelobes of a 
strong source. We see that the arcsecond scale structures have a high contrast with respect to the 
background diffuse emission. We note that it is unlikely that our images can achieve the theoretical 
rms noise level if the extended diffuse emission has dominant contributions to the total flux. 

The small scale emissions may come from compact dense molecular cores, or from local struc- 
tures with high temperature or high CS abundance owing to the interactions with the (proto-) stellar 
activities and shocks. The compact dense molecular cores are potentially the future or current sites 
of star formation. In Figure [71 we see that two significant local maxima of CS (1-0) emission 
are spatially associated with water maser sources (water maser E, and the cluster of water maser 
sources around the UC Hii region A; Hofner & Churchwell 1996). Figure [8] shows the PV diagrams 
of three cuts through the water maser sources N, E, and W. From this figure we see detections of 
bright components near (<1") the water maser sources E and W. The water maser source E may 
be associated with the local maxima of CS (1-0) emission at -7.8 kms~^, or -5.4 kms~^. We do not 
find bright local CS (1-0) emission peak around the water maser source N, which may be due to the 
fact that its parent molecular core is already dispersed by the outfiows. Recalling the small high 
opacity clumps detected in the NH3 (3,3) hyperfine inversion line absorption (Sollins and Ho 2005), 
it is not surprising that more of them are detected in emission line. However, our current data 
are insufficient to robustly distinguish the dense cores from the local shock and outfiow signatures. 
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which leads to the difficulties in a quantitative analysis of the core mass function. In the other 
paper on high velocity molecular outflows in this region, we will discuss the feedback from those 
local sites of star formation and the implications of clustering of stars (Liu, Ho & Zhang 2010). 

3.2.2. Dynamical Motions 

We show the mean velocity map (moment one) and the velocity dispersion map (moment two) 
of the CS (1-0) emission in Figure EJ From Figure El we see that the entire system has a bimodal 
velocity profile. In addition to the general rotation with velocity gradient from the southeast to 
the northwest, we see blueshifted gas in the northeast and redshifted gas in the southwest. These 
additional blueshifted and redshifted gas cover a larger velocity range than the general rotational 
motion, and the orientation resembles a bipolar expansion with a certain inclination. Geometrically, 
the upside-down V shaped cavity-wall signature in the southwest can also be explained by the 
bipolar expansion. We note that our 3(T detection limit is about 18 K, which is only sensitive to 
structures which are sufficiently warm and have high optical depth. The moment one and moment 
two maps thus can be biased due to the cutoff of the fainter high velocity line wings, and over 
interpretations should be omitted. 

To understand the dynamics of the cavity-wall, we make two PV cuts through the southern 
part of the biconical cavity (centered at cutl: R.A. = 18''10"^28.51* , Decl = -19°55'54.3", and 
cut2: R.A. = 18''10™28.48* , Decl = -19°56'1.0") with a position angle 140° (consistent with the 
rotational plane of the central compact hot core PA=140°ib5°; Liu et al. 2010); and we make one 
PV cut through the UC Hii regions B, C, and the northern part of the biconical cavity (centered at 
cutn: R.A. = 18''10'"27.45' , Decl = -19°55'48.8", with position angle 70°). Features we identify 
from these three PV cuts are described as follows. Owing to the complexity of the observed region 
and the limited signal-to-noise ratio of the observations, there could be other dynamical features 
not revealed in the current maps. 

cutn From the PV diagram of the cutn (Figure [TOl) . we identify the expansional signature in the 
northern part of the biconical cavity. Around the UC Hii region B, we marginally see the arc-shaped 
signatures which have broad velocity width, and may also be explained as expansional signatures 
around the UC Hll region B. The expansional velocities in these two regions are about 3 kms~^. 
In these two regions, the expansional signatures are confused with the global rotation/contraction 
of the dense gas, and therefore the size and the expansional velocities are not robustly constrained. 
The detection of the expansional signature in the UC Hii region C is marginal, with a suggested 
expansional velocity of 3.8 kms~^. 

cutl and cut2 Both the cutl and cut2 are located south to the dense flattened structure; the 
center of cut2 is south of the center of cutl. We argue that PV cuts in the position angle of 140° 
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is less confused by the dynamics of the fast rotating (~3-4 kms~^ rotational velocity) dense gas in 
the mid-plane of the envelope, and provides the clean information of the dynamics of the molecular 
cavity-wall. 

The top and bottom panels of Figure \TT\ show the PV diagram of cutl and cut2, respectively. 
The CS (1-0) emission in both panels clearly show arc-shaped expansional signatures centered 
at v;sr~-3 kms~^, with an averaged brightness temperature of 20-30 K, and the peak brightness 
temperature of 45 K. Supposedly CS (1-0) is optically thick at the peak, then the brightness 
temperature represents the excitation temperature of the molecular gas. The maximum excita- 
tion temperature of 45 K in the projected area of the biconical cavity explains why the CH3OH 
transitions are difficult to be excited and have low optical depth. The projected linear scale of 
the expansional signature is about 10" (0.3 pc), which is consistent with the scale of the visually 
identified biconical features in the CS (1-0) maps (Figures [5l E]), and the CH3OH maps (Figures El 
HI). 

The arc-shaped signatures in the PV diagrams have a different physical size scale and geometry 
from the typical outflow signatures (Liu, Ho, & Zhang 2010). The arc-shaped feature shown by the 
cutl may be explained as the redshifted expansional motion of the cavity-wall. Since the center 
of cutl is closer to the free- free continuum peak, its blueshifted counter part may be affected by 
the free-free continuum subtraction, and therefore we cannot obtain significant detections; or, it 
can be explained by the non-uniformity of excitation temperature and density distribution (the la 
detection limit is 5.9 K). The arc-shaped feature shown by cut2 can be explained as the blueshifted 
expansional motion. The lack of the redshifted counter part of the expansional motion can either 
be explained by the asymmetry of molecular gas density, or asymmetrical stellar ionization. The 
line-of-sight expansional velocity of both features is about 2.5-5 kms~^. This interpretation can be 
constrained by observing the lower density tracers (e.g. ^^CO (2-1) or C^^O (2-1) ), to sample the 
lower density/temperature part of the cavity wall. The emissions of those lower density tracers are 
extended, and the single dish observations will be needed in order to robustly interpret the data. 



4. Discussions 

4.1. The Overall Picture of the O Type Cluster Forming Region 

We summarize our series of high resolution observations on GlO.6-0.4 by a schematic model, 
which is shown in Figure [T2l The observational results enable the discussions of the importance 
of the radiation and the pressure force of the ionized gas. In Sections 14.21 and 14. 3^ we provide 
the order-of-magnitude estimates of the radiative pressure and the pressure forces of the ionized 
gas and the stellar wind. We discuss their implications for the dynamical evolution of the system 
in Section 14.41 We are looking forward to comparing our results with numerical hydro dynamical 
simulations in the future. 



- 13 - 



The dominant molecular structure in this model is a 0.5 pc scale massive envelope. The 
massive envelope rotates and contracts globally. A flattened overdensity formed in the mid-plane 
of the massive envelope (this is supported by the CS (1-0) observations, and also the ^'^CS (5- 
4) observations in Liu et al. 2010). Meanwhile, the molecular gas also contracts locally, and 
forms intermediate-mass or B-type massive (proto-)stars, and some low-mass (proto-) stars. The 
protostellar objects eject molecular outflows with high momentum and high energy into the ambient 
environment, which may affect the global contraction of the massive envelope (Li & Nakamura 2006; 
Nakamura & Li 2007; Carroll et al. 2009; Wang et al. 2010). Those molecular outflows can be 
revealed as highly blueshifted and redshifted ^^CO emissions (|u|>20 kms~^) if the protostellar 
objects are not located in the mid-plane (Liu, Ho, & Zhang 2010). In the mid-plane, the outflows 
are impeded and decelerated by the high density ambient gas to have lower terminal velocities, 
which have to be diagnosed from the molecular outflow/shock tracers (Liu, Ho, & Zhang 2010). 

The regulated global contraction may lead to the formation of the 0-type cluster in the center 
(Li & Nakamura 2006), which creates an UC Hii region. The ionized gas and the ionizing photons 
leak out from the Hii region mainly in the bipolar direction, either because of a configuration 
where the rotationally flattened system has a low gas density in the bipolar region, or because of a 
strong protostellar MHD wind which created the biconical cavities in that region. The outflowing 
molecular/ionized gas may flush through the ambient molecular gas on a timescale comparable 
to or much shorter than the global dynamical timescal^, and leave behind the elongated relics 
of interactions. This may explain the detected large fllamentary structures by CS (1-0) in the 
northeast (Section I3.2.ip . Due to the short expansional timescale, the filamentary structures were 
not yet relaxed by the global rotational/infall motion. The pressure of the ionized gas can also 
push the cavity-wall to a quick dynamical expansion (Section I3.2.2p . 

However, in the rotationally flattened system, most of the molecular mass are concentrated in 
the mid-plane. Owing to the self-gravitational or the Rayleigh-Taylor instabilities, the molecular 
gas becomes clumpy (Sollins & Ho 2005), and additionally reduces the effective solid angle as 
presented to the central stars. In such cases, the momentum impulse exerted by the ionized gas 
would be much smaller because of this smaller solid angle. The inertia of the molecular accretion 
flow in the mid-plane can easily overcome the impulse exerted by the stellar radiative pressure 
(also see Section S]) and the pressure of the ionized gas. This is consistent with the lack of a clear 
expansional motion in the P V diagram of the molecular lines (Figure [3]) . When the molecular gas 
gets sufficiently close to the central OB cluster, it is heated by the stellar radiation to a significantly 
higher temperature, and is further ionized to a temperature of 10^ K (also see discussions in Keto 
2002; Keto & Wood 2006). This overall scenario allows the 0-type stars to continue accreting after 
the nuclear burning begins. 



''The lower limit of the global dynamical timescale can be estimated by the free-fall timescale of ~10^ years. 
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4.2. The Role of the Biconical Cavity in Enhancing the Molecular Accretion 

The biconical cavity provides a low opacity channel for the photons to leak out of the envelope. 
This mechanism can significantly reduce the radiative pressure on the molecular gas, and can 
therefore enhance the molecular accretion flow. 

To quantitatively examine this effect, we first compare the relative importance of radiative 
pressure and the gravitational force in a system without the biconical cavity. Following the deriva- 
tions in Jijina &: Adams (1996), we assume the radiation field is nearly isotropic; and the stellar 
mass predominantly contributes to the gravitational force in the relevant scale. We assume the 
temperature distribution follows the power law r(r) = ro(r/ro)-^/2 (Keto, Ho, & Haschick 1987), 
and quote the resulting effective potential from Jijina &: Adams (1996): 



GM 



eff 



r 

where M is the mass of the embedded stellar cluster. The parameter a is defined by 

_ Lkp{To) 

where L is the bolometric luminosity of the embedded stellar cluster, kp{To) is the Planck mean 
opacity at the fiducial temperature Tq, and rg is the radius where T = Tq. The first term in equation 
[T] represents the radiative pressure, and the second term represents the gravitational force. From 
equation [U we see that when a^/r ~ 1, the radiative pressure is comparably important with the 
gravity. The radiative pressure becomes dominant at smaller radii. In GlO.6-0.4, the molecular 



hot toroid with temperature about 300 K is detected inward of the 0.05 pc radiu^j. To match the 
measured molecular gas temperature, we adopl^l Tq = T^; ~ 2300 K (k(To)~30 cm^g~^) and ro = 
8.5-10~^ pc. Given the GlO.6-0.4 bolometric luminosity L ~10^ L0 and stellar mass M ~200 M0, 
we obtain ~0.07 pc (Ho & Haschick 1981; Keto 2002; Sollins & Ho 2005; Keto & Wood 2006). 
In such a case, the radiative pressure can significantly affect the dynamics at the scale of the hot 
toroid, or even reverse an inflow. 

However, based on the numerical radiation transfer calculations of a set of systems with dif- 
ferent geometries, Krumholz, McKee and Klein (2005) suggested that the presence of the biconical 
outflow cavity allows the photons to leak out, and therefore potentially reduces the radiative pres- 
sure force by a factor of ~10. This effectively reduce L by a factor of 10, and therefore reduce 



''The averaged temperature might be overestimated. A lower value of 87 K is reported by Beltran et al. (2011) 
by observing optically thinner lines. We adopt the previous higher value to provide an upper limit of the radiation. 

^We express the temperature power law with a parameter Ta, which is the dust sublimation temperature. If the 
temperature is higher than Td, the dust is sublimated and the radiative pressure is reduced. The temperature power 
law of T(r) = To{r/ro)~^^^ potentially overestimates the temperature at the 0.3 pc radius by a factor of 2 since the 
temperature decay faster with radius in regions with lower opacity, and has to be fitted with a more negative power 
law index. However, our discussions about the radiative pressure focus on the inner region which has high opacity. 
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by a factor of 100. Assuming this conclusion can be generally applied, with the presence of the 
biconical cavity, the value of o? in the case of GlO.6-0.4 might be effectively reduced to 7-10~^ 
pc (140 AU), which is much smaller than the radius of the UC Hii region in GlO.6-0.4 (~0.03 
pc). Qualitatively this estimate is consistent with the case in Krumholz, McKee & Klein (2005), 
suggesting that without considering the ionization, the radiative pressure will only be important 
on the scale of a few hundred AU. We therefore suggest that in GlO.6-0.4, the radiative pressure 
is unimportant for the dynamics of the molecular accretion flow. We note that we may overesti- 
mate the reduction factor of the radiative pressure owing to the non-self-consistent treatment of 
the temperature distribution. However, since o? scales as the square of the reduction factor, it can 
easily become smaller than the radius of the Hii region with a small reduction factor. The presence 
of the biconical cavity can also lead to a temperature gradient steeper than our assumption, and 
may makes o? shrink to an even smaller radius. 

The dusty disk around the stars may redirect the radiation (Yorke & Sonnhalter 2002; Kuiper 
et al. 2010), and potentially play the role of reducing the radiative pressure on the accretion flow. 
However, it is still uncertain whether the disks can stably exist in UC Hii regions, where multiple 
0-type stars are embedded. This has to be examined in future observations. 

4.3. The Ionized Gas Pressure and the Stellar Wind Feedback 

The large molecular mass in the contracting massive envelope contains enormous inertia. Prom 
the PV diagrams (Figure [3]), we find that the rotational velocity at the 0.25 pc radius is about 4-5 
kms~^. This rotational motion is gravitationally bound by an enclosed mass of 1000-1600 Mq. 
While the embedded OB stars contribute ~200 Mq, the rest of the binding mass is dominantly 
contributed by the molecular gas. For a spherical geometry, this molecular mass corresponds to a 
mean molecular density nli^ of (3-5) -10^ cm""^. The momentum flux in the molecular gas flow can 
be estimated by pb-nu^-v^ ^ where v is the infall velocity of the molecular gas, which can be larger 
than 1 kms"^ (Ho & Haschick 1986; Keto, Ho & Haschick 1987, 1988; Keto 1990; Klaassen & Wilson 
2008), and [i is the mean molecular weight. In the embedded UC Hii region, the momentum flux 
in the ionized gas can be estimated by fii-rii-vf, where fii is the ion mass, nj is the ion density, 
and vf is the thermal velocity of the ions, which has the order of magnitude of 10 kms~^. We 
assume that the ion density nj is equal to the electron density ne, which is constrained to be 10^- 
10^ in previous observations (Ho and Haschick 1981; for a more sophisticated model fitting see 
Keto, Zhang, & Kurtz 2008). Assuming //~2/Xi, we find that the momentum flux in the molecular 
gas flow (6-10-10^-^i M0kms~^cm~^s~^) and the momentum flux in the ionized gas (l-10-10^-/ii 
M0kms~^cm~^s~^) have the same order of magnitudes. This suggests that the two components are 
in an approximate dynamical equilibrium. In the real case, the mass is anisotropically distributed. 
The dynamics of the accretion flow in the denser region will not be significantly affected by the 
pressure of the ionized gas, while the lower density region can undergo a dynamical expansion 
powered by the pressure of the ionized gas. 
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We estimate the momentum budget of the ionized gas and the stellar wind in this section. In 
the next section, we will suggest a simple geometric model for the density distribution, and estimate 
the velocity of the dynamical expansion, to compare with the observations. The embedded massive 
stars in UC Hii regions B and C have similar spectral type (Ho & Haschick 1981) with the embedded 
OB stars in UC Hii region A. Some of the energetic relations estimated in UC Hli region A can 
also be applied in UC Hii region B and C. 

Expansional Time Scale The expansional signatures of the biconical cavity have a character- 
istic radius rc of about 5". Assuming the fastest expanding front of the cavity- wall uniformly 
expands with the terminal velocity Ve of 5 kms~^, the total expanding time tg is about 3T0^ years, 
which is much shorter than the global dynamical timescale (>10^ years). The timescale estimation 
constrains the total momentum budget from the stellar wind and the ionized gas. 

The Stellar Wind Feedback The UC Hii region A contains a few 06-09 stars with a total 
stellar mass of ~200 M©. The stellar wind from each of these massive stars has a mass loss 
rate of the order of magnitude of 10~^ MQyr"^, and terminal velocity of ~2000 kms~^ (Tout et 
al. 1996). Given the number of the embedded O stars (~4; Ho & Haschick 1981), the wind 
momentum feedback rate of the massive cluster is A^*-2-10~^ MQkms~^yr~"'^. Assuming the 
embedded stellar cluster starts to feed back the stellar wind by the beginning of the expansion of 
the biconical cavity, the total wind momentum budget I^j is N^:-60 MQkms"^ isotropically spread 
in the entire Att solid angle. The exerted pressure force has a radial dependence. 

The Ionized Gas Pressure We adopt the characteristic radius r of 0.15 pc, and adopt the 
electron density rig = 10^-10^ cm~^ according to the measurements in Ho & Haschick (1981), to 
estimate the thermal pressure force of the ionized gas. Assuming the ion density equals to the 
electron density, and has the thermal velocity of ~10 kms^^, the momentum feedback rate fi can 
be estimated by Q-r'^ -{p-v -v), where is the characteristic solid angle. As an order of magnitude 
estimation, we adopt 0=47r, which leads to /j~6T0~^-6T0~^ MQkms~^yr~^. The total ionized 
gas momentum feedback = f^" fi{r{t))dt depends on the expansion history, and has the order of 
mag nitude of 10^-10^ Mokms^^ 

The ionized gas feedback differs from the stellar wind feedback in the sense that it is a pressure 
effect (see also Figure [T2]) . which does not exert force in the radial direction, but in the direction 
perpendicular to the cavity-wall. In addition, the exerted pressure force depends only on the 
temperature and density of the ion, which is determined by the ionization, the recombination, and 
the pressure balance. 
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4.4. The Dynamical Process 

From the observations of the 1.3 mm continuum emission, and the measurements of the ro- 
tational velocity of the molecular gas, we constrained the molecular mass inward of the 0.15 pc 
radius around the UC Hii region A to have the order of magnitude of 400 Mq . We assume a simple 
geometry in this region, and perform order-of-magnitude estimates of the energetic relation. 

4-4- 1- The Geometrical Model 

Assume the molecular gas has an initial mass density distribution po before the creation and the 
dynamical expansion of the biconical cavity. Defining 77 to be the angular separation from the plane 
of rotation, if the system is approximately axisymmetric, po can be represented by po = Po{i"i v)- If 
Po{r) can be represented as a polynomial of cos(ry), for a rotationally flattened system, we expect 
the leading order dependence to be cos(r/). This suggests that about 70% of molecular mass (280 
Mq) are concentrated in the region with |r/| <45°. 

After the creation and the dynamical expansion of the biconical cavity, the molecular mass 
initially distributed in the region \rj\ >45° (120 Mq) are accumulated on the geometrically thin 
cavity-wall, which has an opening angle of 90°. The assumption of the mass accumulation on the 
cavity-wall is consistent with the fast expansion (see Section 14.31 for the estimate of timescale) . 
The molecular gas initially distributed in the region \r]\ <45° forms a dense flattened structure 
with density distribution p. We argue that the opening angle of 90° is a reasonable value while 
comparing with the geometrical picture revealed by the CH3OH images (Figure [21 H]) . The cavity- 
wall has two components, W_l and Wy, of which the surface area are perpendicular and parallel 
to the radial direction, respectively. Given the opening angle, the ratio of the surface area of these 
two components is approximately 1; we assume the ratio of the mass accumulated on these two 
components is also 1. 

4-4-^- The Dynamical Expansion of Cavity Wall 

W|| This component has a cone shape, with the total mass of ~60 Mq. For it to expand at 
2.5-5 kms^^, it requires a momentum of 150-300 MQkms"^. From the estimations in the previous 
section, we see that this required momentum has the same order of magnitude of the feedback from 
the ionized gas pressure. The gravitational force acts in the radial direction and does not retard 
the expansion in the 6 direction. We expect the stellar wind to have weak effects on the dynamics 
of W|| since the wind is parallel to the cavity-wall. 

This part is directly driven by the stellar radiation and the stellar wind, which are competing 
with gravity in the radial direction. At the radius of 0.15 pc, the stellar wind and the ionized gas 
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pressure can have comparable effect on W_l. Assuming the mass is 60 Mq, the gravitational force 
has the order of magnitude of 2-10"'^ MQkms~^yr~^, which is also comparable to the force of the 
stellar wind and the ionized gas. 

Observationally this part of cavity-wall is not as prominently detected as Wy. Since the 
embedded stellar cluster in the UC Hii region A has a high ionizing photon emission rate {S = 10^^ 
s~ ; Ho & Haschick 1981), which is able to ionize on the order of 10^ Mq of gas in the biconical 
region. Even if the recombination rate marginally balances the ionization rate, the stellar ionization 
can still explain why the is not clearly detected. Some local high density clumps which have 
higher recombination rate may be self-shielded from the stellar ionization and have long lifetime. 
Those high density clumps occupy small solid angles, and are less accelerated by the stellar wind 
and the pressure of the ionized gas. Some of those clumps can still fall toward the OB cluster 
owing to the gravitational attraction. The clumpy 3.6 cm free-free continuum emissions in the 
UC Hii region A, especially, northeastern to the emission peak, may be explained by those clumps 
externally ionized by the central OB cluster. 

We note that we may overestimate the dynamical expansion timescale of the biconical cavity 
since the cavity can start expansion from a finite size owing to the initial ionization, which leads to 
the overestimation of the momentum feedback. However, those ionized gases are not gravitationally 
bound and are free to leak out. Therefore, the accumulated mass on the cavity-walls and Wj_ 
is reduced if the initial ionization is efficient, and requires smaller momentum feedback to reach the 
final state velocity. Similar arguments are valid if the biconical cavity is initially created by the 
massive bipolar outflow (an example see G240.31+0.07: Qiu et al. 2009) and the molecular gas is 
evacuated. 

The Dense Flattened Structure This structure has the mass of ~280 Mq, and continues from 
the outer of ~0.15 pc radius to the inner ~0.03 pc (Liu et al. 2010). The total gravitational force 
on this structure depends on the embedded mass. A lower limit of 9.6T0~^ MQkms"^yr"^ can be 
given assuming that the majority of mass is distributed at the 0.15 pc radius. This lower limit 
is already on the same order of magnitude as the feedback from the stellar wind and the ionized 
gas pressure; the actual value of the gravitational force should be much larger. Depending on the 
effective solid angle it occupies, the stellar wind and the ionized gas pressure may have negligible 
effect on the radial motion of the flattened structure. With the biconical cavity structure which 
leads to the photon leakage, the radiative pressure is also negligible. We therefore suggest that the 
dynamics in this region is dominated by gravity and rotation. 

UC HII region B and C The stellar wind and ionized gas pressure feedback in UC Hii region 
B and C have the same order of magnitude as those in the UC Hii region A. However, these two 
regions have lower initial gas density, which leads to a smaller accumulated mass on the expansional 
shell. The observed expansional velocities in these two regions are comparable to but smaller than 
the thermal velocity of the ionized gas, suggesting that the ionized gas pressure is an important 



- 19 - 



driving source of the dynamical expansion. From the estimations in the previous section, the stellar 
wind feedback can be comparably important, depending on the detailed geometry of these system. 

5. Summary 

We present high resolution observations of molecular lines and free-free continuum for the UC 
Hii region G10.6~0.4. The resolved projected distributions of the molecular gas and the ionized gas 
suggest an overall picture consisting of an extended (~0.5 pc) envelope, a single compact (~0.1 pc) 
hot rotating toroid, and a biconical molecular cavity filled with ionized gas. This overall geometry 
resembles the standard envelope-disk and protostellar outflow model for the low-mass star forming 
region. With the presence of the biconical cavity, we suggest that the radiative pressure can be 
significantly reduced. In the plane of rotation, at the scale larger than 1" (0.03 pc), we see that the 
rotational motion of the dense gas is not severely disturbed, which consistently suggests that the 
radiation is not yet important for the molecular accretion flow. The stellar radiation may increase 
and play a more important role as the stellar mass is increased. 

From the observations of CS (1-0), we suggest that the biconical cavity around the UC Hii 
region A is undergoing an expansional motion, with velocity of 2.5-5 kms~^. We perform simple 
order of magnitude estimates, and suggest that the feedback from the ionized gas pressure can 
account for the required momentum of this dynamical expansion. The expansional signatures 
are also detected in the UC Hll region B and C, which are driven by the ionized gas pressure 
and the stellar wind. The expansional motions of the UC Hli regions significantly disturb the 
local dynamics of the molecular gas, inject the energy, and may induce the non-uniformity in the 
molecular accretion flow, which is an important feedback mechanism in the massive molecular 
clump. 

Facilities: SMA, VLA/EVLA 
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Transition 


Frequency (GHz) 


E„p/k (K) 


Eins. A. (s^i) 


Note 


NH3 (3,3) main 


23.870129 


124.5 


2.56-10- 


7 




CS (1-0) 


48.990055 


2.35 


1.75-10" 


6 




CH3OH 5(0,5)-4(0,4) E 


241.700219 


47.68 


6.04-10- 


5 




CH3OH 5(0,5)-4(0,4) A+ 


241.791431 


34.65 


6.05-10" 


5 


R 


CH3OH 5(-2,4)-4(-2,3) E 


241.904152 


60.38 


5.09-10" 


5 


G 


CH3OH 5 (2,3)-4(2,2) E 


241.904645 


57.27 


5.03-10" 


5 




CH3OH 5(-3,3)-4(-3,2) E 


241.852352 


96.93 


3.89-10" 


5 


B 



Table 1:: Table of the selected molecular transitions. The quantum number of the transitions are 
listed in the first column. Their frequencies and upper-level energy are listed in the second and the 
third column. The Einstein A~Coefficient of each transition is listed in the fourth column. The 
last column notes the corresponding color of the CH3OH lines in the RGB images (Figure [2l 



Transition Synthesized beam RMS noise (Jy/beam) Instruments uv sampling range (kX) Observed Date 



NH3 (3,3) main 


1" 


8x1". 2 


0.008 


VLA/EVLA C-array 


2.2-270 


2009.07.27 


CS (1-0) 


1" 


5x1". 1 


0.022 


VLA/EVLA DnC-array 


4-245 


2009.09.27 


CH3OH 5(0,5)-4(0,4) E 


1" 


5x1". 3 


0.06 


SMA compact -t-very extended 


6-393 


2009.06.10/2009.07.12 


CH3OH 5(0,5)-4(0,4) A-l- 


1" 


5x1". 3 


0.06 


SMA compact -f-vcry extended 


6-393 


2009.06.10/2009.07.12 


CH3OH 5(-2,4)-4(-2,3) E 


1" 


5x1". 3 


0.06 


SMA compact+very extended 


6-393 


2009.06.10/2009.07.12 


CH3OH 5 (2,3)-4(2,2) E 


1" 


5x1". 3 


0.06 


SMA compact -|-very extended 


6-393 




CH3OH 5(-3,3)-4(-3,2) E 


1" 


5x1". 3 


0.06 


SMA compact+very extended 


6-393 


2009.06.10/2009.07.12 



Table 2:: Table of the instrumental parameters of the selected molecular transitions. The NH3 
data has the velocity resolution of 1.2 kms~^; the CS (1-0) data has the velocity resolution of 0.6 
kms~^; and the CH3OH data has the velocity resolution of 0.5 kms~^. The rms noises shown in 
the third column are measured from the channel images with corresponding velocity widths. 



Observation Flux Cal. Passband Gal Gain Cal Total Bandwidth Frequency Resolution 



VLA/EVLA G-array (NH3) 


1331+105 


3c273 


1820-254 




6.25 MHz 


97.656 kHz 


VLA/EVLA DnC-array (CS) 


1331+105 


0319+415 


1733-130/1820-254 




6.25 MHz 


97.656 kHz 


SMA compact-array 


Uranus 


3c273 


1733-130/1911-201 


2 GHz 


in each sideband 


406 kHz 


SMA very-extended-array 


mwc349 


3c273 


1733-130/1911-201 


2 GHz 


in each sideband 


406 kHz 



Table 3:: The observational settings of the selected molecular transitions. 
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Fig. 1. — : The velocity integrated maps of the CH3OH J=5 transitions (Left: CH3OH 5(0,5)- 
4(0,4) A+; Middle: the blended CH3OH 5(-2,4)-4(-2,3) E and CH3OH 5 (2,3)~4(2,2) E; Right: 
CH3OH 5(-3,3)-4(-3,2) E). Color-bars have the unit of Jy/beam-kms~^. Three 1.3 cm free-free con- 
tinuum peaks (R.A.=18''10™28^683 Decl=-19°55'49".07; R.A.=18''10'"28^215 Decl=-19°55'44".07 
; R.A.=18''10'"27^435 Decl=-19°55'44".67) are marked by red stars. The synthesized beams are 
shown in the bottom left corner of the panels. 
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Fig. 2.—: Left: The RGB image of the CH3OH J=5 transitions (R: CH3OH 5(0,5)-4(0,4) A+; 
G: the blended CH3OH 5(-2,4)-4(-2,3) E and CH3OH 5 (2,3)-4(2,2) E; B: CH3OH 5(-3,3)-4(-3,2) 
E). Right: The RGB image of the CH3OH J=5 transitions overlaid with the NH3 (3,3) main 
hyperfine inversion emission (contour). Solid contom's start from 10% of the emission peak with 
10% intervals; dashed contom's start from -200% of the emission peak with -200% intervals. Note 
the primary beams of the SMA (~l') and the VLA (~2') observations are much larger than the 
angular size scale of the detected structures. 
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Velocity (km/s) 



Fig. 3. — : The PV diagrams of the CH3OH J=5 transitions (grey scale and black contours) cut 
in the plane of rotation, which is centered at the coordinates of R.A. = IS'* 10™ 28.64'' and Decl 
= -19°55'49.22" with position angle pa = 140° From left to right is the CH3OH 5(0,5)-4(0,4) 
A+ transition, the CH3OH 5(0,5)-4(0,4) E transition, the blended CH3OH 5(-2,4)-4(-2,3) E and 
CH3OH 5 (2,3)-4(2,2) E transitions, and the CH3OH 5(-3,3)-4(-3,2) E. The black contour intervals 
and the first contour level are 0.18 Jy/beam (1.9 K) for each panel. In the right most panel, we 
also plot the NH3 (3,3) main hyperfine emission line PV diagram in blue contours, start from 0.01 
Jy/beam (11 K) with 0.01 Jy/beam intervals; and plot the NH3 (3,3) main hyperfine absorption 
line by one dashed green contour at the level of -0.5 Jy/beam. We plot the 0.01, 0.04, and -0.5 
Jy/beam contour levels of the NH3 (3,3) main hyperfine line in the left most panel for comparison. 
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Fig. 4. — : The RGB image of the CH3OH J=5 transitions overlayed with the 3.6 cm continuum 
image (contour). Contours show the levels of [1, 2, 4, 8, 16, 32, 64, 96]% of the peak intensity of 
0.1 Jy/beam. 
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Fig. 5. — : Grey Scale: The moment map of the CS (1-0) emission. The synthesized beam of 
this map is 1.4" x 1.2", with position angle 77°. Contour: The 1.3 mm continuum image. This 
data is taken by SMA subcompact+compact+very extended array with 0". 79x0". 58 resolution. 
Contours show the levels of [1, 2, 4, 8, 16, 32, 64, 96]% of the peak brightness temperature 49 
K. A proposed plane of rotation is indicated by a red dashed line. The red dashed line passes 
through two free-free continuum peaks, and has a position angle of 135°, which is consistent with 
the position angle of the plane of rotation of the central 0.1 pc scale hot core (140°ib5°, Liu et al. 
2010). 
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Fig. 6. — : Grey Scale: The moment map of the CS (1-0) emission. The synthesized beam of 
this map is 1.4" x 1.2", with position angle 77°. Contour: The 3.6 cm free-free continuum image; 
resolution 0".5x0".4. Contours show the levels of [1, 2, 4, 8, 16, 32, 64, 96]% of the peak intensity 
of 0.1 Jy/beam. 
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Fig. 7. — : The example channel map at Yisr= -7.8 kms ^. Contours start from -0.06 Jy/beam 
with 0.06 Jy/beam intervals. The synthesized beam is shown in the bottom left corner of 
this image. Three 1.3 cm free-free continuum peaks (R.A.=18''10™28^683 Decl=-19°55'49".07; 
R.A.=18''10""28^215 Decl=-19°55'44".07 ; R.A.=18^10'"27^435 Decl=-19°55'44".67) are marked 
by red stars. Dot symbols mark the water maser detections (Hofner and Churchwell 1996). Three 
water maser sources which are not spatially associated with the UC Hii region A are labeled by N, 
W, and E, respectively. The relative positional accuracy of the maser data is about 0.1", which is 
smaller than the size of the Dots. 
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Fig. 8. — : From left to right, the position-velocity (PV) diagrams at the locations of water maser 
sources N, E, and W. Contours start from -0.06 Jy/beam with 0.06 Jy/beam intervals. The syn- 
thesized beam is shown in the bottom left corner of this image. Positive angular offset is defined 
in the east. The PV cuts center at the locations of the water maser sources, and have the position 
angle PA=90°. Color-bars have the unit of Jy/beam. The deep dark region in the right most panel 
is the absorption line against the bright background continuum emission. 
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RA (J2000) RA (J2000) 

Fig. 9. — : The mean velocity map (left) and the velocity dispersion map (right) of the 
CS (1-0) emission. Color-bars have the unit of kms~^. Three 1.3 cm free-free contin- 
uum peaks (R.A.=18''10™28^683 Decl=-19''55'49".07; R.A.=18''10'"28^215 Decl=-19''55'44".07 
; R.A.=18^10'"27^435 Decl=-19°55'44".67) are marked by blue stars. Dot symbols mark the water 
maser detections (Hofner and Churchwell 1996). Three water maser sources which are not spa- 
tially associated with the UC Hii region A are labeled by N, W, and E, respectively. The relative 
positional accuracy of the maser data is about 0.1", which is smaller than the size of the Dots. 
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Fig. 10. — : The PV diagram cuts through the UC Hii region B and C. The center of the PV cut is: 
R.A. = 18''10""28.51* , Decl = -19°55'54.3"; the position angle of the cut is 70°. Positive angular 
offset is defined in northeast. Solid contours in both panels start from the brightness temperature 
of IIK (~2(7), with IIK intervals; negative contours start from the brightness temperature of - 
IIK, with -IIK intervals. Blue circles from top to bottom are indicative of expansional signatures 
with centroid v;jjj,=-4, -1 and -2 kms~^, respectively; and the expansional velocity 3.1, 2.7, and 3.8 
kms~^, respectively. 
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Fig. 11. — : Top: The PV diagram cuts through the southwest lobe (0" angular offset corresponds 
to R.A. = 18^10'"28.51' , Decl = -19''55'54.3") with position angle 140°. Positive angular offset 
is defined in southeast. Bottom: The PV diagram cuts through the southwest lobe (0" angular 
offset corresponds to R.A. = 18^10'"28.48'' , Decl = -19°56'1.0") with position angle 140°. Solid 
contours in both panels start from the brightness temperature of IIK, with IIK intervals; negative 
contours start from the brightness temperature of -IIK, with -IIK intervals. Blue circles are 
indicative of an expansional signature with centroid v;5j.=-3kms~^, and the expansional velocity 
2.5 and 5.0 kms~^, respectively. 
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Fig. 12. — : A schematic view of the 0-type cluster forming massive clump. The yellow star symbols 
indicates the (proto-)stellar objects. A compact 0-type cluster is embedded in the center of the 
system, and is the dominant source of the stellar radiation. The Red and blue arrows indicate 
molecular outflows. Longer arrows represent outflows with larger terminal velocities; the outflow 
momenta are proportional to the area of the arrows. We note that the molecular structures in this 
figure are highly clumpy /filamentary. 



